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We reconsider the model of the randomly pinned planar flux array (or 
the two-dimensional XY-model with random fields and no vortices) in 
extensive Monte Carlo simulations. We show that contrary to what has 
been reported recently [Phys. Rev. Lett. 72, 4133 (1994)] the disorder- 
averaged displacement correlation function is different from the pure case 
for larger values of the disorder. Its slope as a function of the logarithm 
of disorder is compatible with newer RG-calculation results. 
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In a recent letter Batrouni and Hwa Jl| reported results of numerical simulations of the 
planar flux array described by the random phase model 

H = - J2( u i - u if ~ cos(2nUi - p { ) , (1) 

(ij) i 

where (ij) denote nearest neighbor sites on a square lattice, U{ a real- valued displacement- 
like field, Pi a random phase uniformly distributed in the interval [0, 2tt], and A the strength 
of the pinning potential. The main conclusion of the investigation in Q was a) that the 
various renormalization group predictions existing up to that time could be ruled out by 
their numerical results, which I agree with, and b) that the disorder averaged correlation 
function C{r) = [((«, — u i+r ) 2 )] wi where (• • •) means the thermodynamic expectation value 
and [• • -] av the disorder average, is indistinguishable from the pure case, i.e. C(r) = C\ = o(r) = 
T/nTclogr for T <T g = k/tc. 

In this comment I would like to point out that the last statement is a consequence of 
the weakness of the disorder they used, namely A = 0.15, and that for stronger disorder (or 
larger length scales) the correlation function C(r) differs significantly from the pure case. 
Furthermore the numerical data I obtain are compatible with the analytic predictions of 
|||3|] and the recent numerical results of a related model Q with infinite disorder. 

I used the usual Monte Carlo algorithm to calculate the static expectation values of the 
model ([!]). By setting k = 2 the glass transition temperature is T g = 2/n « 0.637. The 
system sizes were L=32 and 64 (where 1280 and 256 samples were used, respectively) and 
periodic boundary conditions were imposed, which forces the correlation function in the 
x-(and y)-direction to be symmetric C(r) = C(L — r), r < L/2. By applying a least square 
fit of the data points for C(r) with r < L/4 (where the effect of the periodic boundary 
conditions is still negligible) to the function C(r) = a + b(T) logr (T being the temperature) 
one gets the results shown in figure 1 for different disorder strengths. 

For weak disorder A ~ 0.15 the slope b(T) is indeed indistiguishable from the pure case 
b(T) = T/ktt, as Batrouni and Hwa observed UJ. However, by increasing the disorder 
(A > 0.5) one obtains a slope b(T) that is significantly different from the pure case already 
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in the vicinity of T g . Since systems with larger disorder are hard to equilibrate, only data 
for A < 2 are shown, but the trend seems to be obvious: the estimate of b(T) obtained from 
intermediate length scales increases with increasing disorder strength. 

One should point out that b(T) shown in figure 1 can only be interpreted as the slope 
of C(r) as a function of logr for small length scales. For certain values of T and A one 
observes a significant upward bend in these curves already at length scale that smaller than 
L/4, indicating that the asymptotic behavior is not yet reached. In figure 2 the derivative 
of C(r) with respect to logr, thus yielding a local slope b(T,r), is depicted for A = 0.5 and 
T = 0.4. This plot gives further support for the hypothesis that lim^oo b(T, r) = T 5 /27r, 
however, a quadratic dependency C(r) ~ log 2 r, implying lim^oo b(T, r) = oo, cannot be 
excluded. With this observation in mind the data shown in figure 1 have to be interpreted 
as lower bounds for the asymptotic slope b(T) of C(r) as function of logr. Hence the results 
for the correlation function are compatible with b(T) = Tg/nrc for T < T g , however, also 
a quadratic dependency C(r) ~ log 2 r, implying lim^oo b(T, r) = oo, although unlikely, 
cannot be strictly excluded. 

Concluding I have presented numerical evidence that the disorder averaged spatial cor- 
relation function of model (|I|) is indeed distinct from the pure case A = 0. For weak disorder 
this becomes manifest only on length scales that are not attainable via Monte Carlo simula- 
tions yet, which is the reason why Batrouni and Hwa ]l| were not able to detect it. However, 
my results agree with their conclusion that the various RG prdictions, prior to their work, 
were incorrect. 

I would like to thank T. Nattermann, S. Scheidl, L. H. Tang and J. Kierfeld for many valu- 
able suggestions and discussions. The computations were done on the Parsytec-GCell024 
from the Center of Parallel Computing (ZPR) in Koln and the Intel Paragon System from the 
Supercomputercenter (HLRZ) at the Forschungszentrum Jiilich. This work was performed 
within the SFB 341 Koln-Aachen-Jiilich. 
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FIG. 1. The slope b(T), obtained via a least square fit of the MC-data for r < L/4, as a 
function of the temperature T. From top to bottom one has A = 2.0 (•), 1.0 (o), 0.5 (o) and 0.25 
(□). Data for A = 0.15, as obtained by Batrouni and Hwa [1] fall onto the thick full line, which 
is the analytic result for the pure model: b(T) = T/ktt. The dash-dotted line indicates the glass 
transition temperature T g = k/tt and the dotted line is at b(T) = Tg/nn. 
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FIG. 2. The local slope 6(T,r) = [C(r + 1) - C(r)]/[Zogi(r + 1) - Zoc/(r)] as a function of 
distance. Shown are data for three different equilibration times: t w = 3- 10 4 (□), i w = 10 5 (+) and 
t w = 3 ■ 10 5 (o). Curves for larger equilibration times are identical to the last one meaning that 
equilibration has been achieved. The dotted line is at b = T g / (ttk) = 1/tt 2 . Because of the periodic 
boundary conditions of the finite system of length L one has b(T,L/2) = 0. If the prediction 
C(r) = Tg/ (ttk) log |r| for T <T g holds, b(T, r) has to approach this line for r — ► 00 in the infinite 
system. 
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